function score= numscore(lambda_0,  rho, Deltas, default, probmf, c1, c2,  param, dx)
if param==1;
    lambda_1=lambda_0;
    lambda_1(1)=lambda_0(1)+dx;
f0= indiv_likelihood(lambda_0, rho, Deltas, default, probmf,  c1, c2);
f1= indiv_likelihood(lambda_1,  rho, Deltas, default, probmf,  c1, c2);
score=(f1-f0)./dx;
end 

if param==2;
    lambda_1=lambda_0;
    lambda_1(2)=lambda_0(2) +dx;
f0= indiv_likelihood(lambda_0, rho, Deltas, default, probmf, c1, c2);
f1= indiv_likelihood(lambda_1,  rho, Deltas, default, probmf,  c1, c2); 
score= (f1-f0)./dx;
end

if param==3;
    lambda_1=lambda_0;
    lambda_1(3)=lambda_0(3) +dx;
f0= indiv_likelihood(lambda_0, rho, Deltas, default, probmf, c1, c2);
f1= indiv_likelihood(lambda_1,  rho, Deltas, default, probmf, c1, c2); 
score= (f1-f0)./dx;
end

    